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ABSTRACT 

We report on the latest (2007 Jan) observations of supernova remnant (SNR) 1987A from the 
XMM-Newton mission. Since the 2003 May observations of Haberl et al. (2006), 11 emission lines 
have experienced increases in flux by factors ~ 3 to 10 (6 ± 0.6 on average), with the 775 eV line of 
O viii showing the greatest increase. Overall, we are able to make Gaussian fits to 17 emission lines 
in the RGS spectra and obtain line fluxes; we have observed 6 lines of Fe XVII and Fe xvin previously 
unreported by XMM-Newton. A two-shock model representing plasmas in non-equilibrium ionization 
is fitted to the EPIC-pn spectra, yielding temperatures of ~ 0.4 and ~ 3 keV, as well as elemental 
abundances for N, O, Ne, Mg, Si, S and Fe. 

We demonstrate that the abundance ratio of N and O can be constrained to < 20% accuracy 
(N/O = 1.17 ± 0.20). Within the same confidence interval, the same analysis suggests that the 
C+N+O abundance varies from ~ 1.1 to 1.4 x 10~ 4 , verifying the Chandra finding by Zhekov et 
al. (2006) that the C+N+O abundance is lower by a factor ~ 2 compared to the value obtained in the 
optical/ultraviolet study by Lundqvist & Fransson (1996). Normalizing our obtained abundances by 
the Large Magellanic Cloud (LMC) values of Hughes, Hayashi & Koyama (1998), we find that O, Ne, 
Mg and Fe are under-abundant, while Si and S are over-abundant, consistent with the findings of As- 
chenbach (2007). Such a result has implications for both the single-star and binary accretion/merger 
models for the progenitor of SNR 1987A. In the context of the binary merger scenario proposed by 
Morris & Podsiadlowski (2006, 2007), material forming the inner, equatorial ring was expelled after 
the merger, implying that either our derived Fe abundance is inconsistent with typical LMC values 
or that iron is under-abundant at the site of Sanduleak -69°202. 

Subject headings: circumstellar matter — methods: data analysis — plasmas — shock waves — super- 
novae: individual (SN 1987A) — ISM: supernova remnants — X-rays: individual 
(SN 1987A) 



1. INTRODUCTION 

SNR 1987A is the Rosetta Stone (Allen 1960) of 
Type II supernova remnants, resolved and well-studied 
in multiple wavebands, including the infrared (Bouchct 
et al. 2006; Kjaer et al. 2007), optical/ultraviolet 
(Groningsson et al. 2007; Heng 2007) and radio 
(Gaensler et al. 2007). The physical mechanism par- 
tially powering optical/ultraviolet emission from the re- 
verse shock is the same as the one at work in Balmer- 
dominated supernova remnants (Heng & McCray 2007; 
Heng et al. 2007). The detection of a neutrino burst 
confirmed the core collapse nature of the progenitor 
(Koshiba et al. 1987; Svoboda et al. 1987), though a 
pulsar has yet to be detected (Manchester 2007). A sys- 
tem of three rings may be the result of a binary merger 
between two massive stars about 20,000 years prior to the 
supernova explosion (Morris & Podsiadlowski 2006, 2007; 
hereafter MP0607). Reviews of the multi- wavelength 
studies of SNR 1987A can be found in McCray (1993, 
2005, 2007). 

Mixing of the stellar envelope and core by Rayleigh- 
Taylor instabilities within the progenitor star, Sanduleak 
-69°202, has been invoked to explain the early emergence 
of the 847 keV 7-ray line from SNR 1987A, which was 
predicted by Shibazaki & Ebisuzaki (1988) to reach its 
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peak around 1.1 years after the explosion, if one assumes 
a mixed mass of about 5M Q (Ebisuzaki & Shibazaki 
1988). Instead, Matz et al. (1988) observed the 847 
keV line ~ 6 to 8 months post-explosion, suggesting even 
more extensive mixing of 56 Co than assumed. A similar 
explanation (Ebisuzaki & Shibazaki 1988) was given for 
the early emergence of 16 to 28 keV X-rays (Sunyaev et 
al. 1990; Inoue et al. 1991). The 7-rays originate from 
the radioactive decay of 56 Co, while the X-rays are from 
the Compton degradation of the 7-rays (McCray, Shull 
& Sutherland 1987). 

In the soft X-ray, SNR 1987A was first observed by 
Beuermann, Brandt & Pietsch (1994). Subsequently, 
Hasinger, Aschenbach & Triimper (1996) tracked a 
steady increase of the soft X-ray flux over 4 years with 
ROSAT. Extensive work has since been done by the 
Chandra (Michael et al. 2002; Park et al. 2002, 2004, 
2005 [hereafter P05], 2006, 2007 [hereafter P07]; Zhekov 
et al. 2005, 2006 [hereafter Z06]) and XMM-Newton 
groups (Haberl et al. 2006, hereafter H06; Aschenbach 
2007, hereafter A07). The general picture gleaned from 
these studies is of a bimodal plasma distribution present 
in the region between the forward and reverse shocks 
(Fig. [J). The soft X-rays (~ 0.3 to 0.5 keV) are from the 
decelerated shock front interacting with dense protru- 
sions ("fingers") on the inner, equatorial ring, while the 
"hard" X-rays (~ 2 to 3 keV) are from a fast shock prop- 
agating into more tenuous material. The soft X-rays ap- 
pear to be correlated with optical "hot spots" , believed to 
be emission from the shocked fingers, appearing around 
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the equatorial ring, while the hard X-ray and radio im- 
ages exhibit structures that coincide. Between Days 6000 
and 6200, the soft X-ray light curve experienced an up- 
turn and departure from an exponentially increasing pro- 
file, which P05 interpreted as evidence that the blast 
wave had reached the main body of the dense circum- 
stellar material of the equatorial ring. Future studies on 
the nature of the soft X-ray light curve are relevant to 
the issue of pre-ionization of the supernova ejecta, which 
can potentially extinguish Ha and Lya emission from the 
reverse shock (Smith et al. 2005; Heng et al. 2006). 

In the studies described, relatively little attention has 
been paid to the subject of inferring elemental abun- 
dances from X-ray analyses. Fits to the X-ray spectra 
yield N, O, Ne, Mg, S, Si and Fe abundances; such re- 
sults have been tabulated and studied by Z06 and H06, 
using Chandra and XMM, respectively. Here, we exam- 
ine such an approach using the latest XMM data set of 
SNR 1987A, taken in early 2007. In SJ1 we describe our 
observations and data reduction techniques. Our results 
are presented in In 21 we perform a detailed error 
analysis of individual abundances and their ratios. We 
find that the N/O ratio as well as the individual N and 
O abundances can be constrained to ~ 20% accuracy. 
Our derived elemental abundances for O, Ne, Mg and Fc 
are under-abundant, while Si and S are over-abundant, 
relative to typical values for the Large Magellanic Cloud 
(LMC; Hughes, Hayashi & Koyama 1998). Such a re- 
sult has implications for modeling the progenitor of SNR 
1987A, which we discuss. 

2. OBSERVATIONS & DATA REDUCTION 

The latest XMM-Newton (Jansen, Lumb & Altieri 
2001) observations of SNR 1987A were performed from 
2007 Jan 17 18:23 to Jan 19 01:18 UT. For our spectral 
analysis, we utilize data from the CCD cameras EPIC-pn 
(Striider et al. 2001a,b) and EPIC-MOS (Turner et al. 
2001), and the Reflection Grating Spectrometers (RGS; 
den Herder et al. 2001). Additional details and net expo- 
sure times for the different instruments are summarized 
in Table [TJ 

The data were processed with the XMM-Newton Sci- 
ence Analysis Software (SAS) version 7.1.0. Using ds9, 
we extracted source and background EPIC spectra of 
SNR 1987A from circular regions placed on the source 
and a proximate, point source- free area (Fig. [5]). For the 
EPIC-pn spectra, single-pixel (PATTERN 0) events were 
selected; for EPIC-MOS all valid event patterns (PAT- 
TERN to 12) were used. RGS spectra were produced 
using rgsproc. The spectra were binned to contain a 
minimum of 20 and 30 counts per bin for EPIC and RGS, 
respectively. 

3. RESULTS 

3.1. SPECTRAL FITTING 

Spectral fitting was done using XSPEC (Arnaud 1996) 
version 11.3.2. We fit the reduced spectra using models 
with two components: a low (Ti ow ) and a high (Thigh) 
temperature component. Two models are considered in 
XSPEC: VNEI+VRAYMDND and VPSH0CK+VPSH0CK. In H06, 
the former ("Model A") is seen to give an excellent 
fit ("reduced chi-square" of \r ~ 1-1) to the EPIC-pn 
data. However, the VRAYM0ND component of the model, 
which describes a plasma in ionization equilibrium, is 



at odds with the belief that both plasmas are in non- 
equilibrium ionization (NEI; Z06); we note that Park 
et al. (2004, 2006) found the low-temperature plasma 
component to be in a highly-advanced ionization state 
(with ionization ages ~ 10 12 to 10 13 cm -3 s) and can be 
described by a thermal plasma in collisional ionization 
equilibrium. The VNEI model is a somewhat unsatisfac- 
tory representation of the physical situation in the post- 
shock plasma of SNR 1987A, because it considers only 
one value of the ionization age/time, r, for the entire 
plasma. The VPSH0CK model generalizes the VNEI one, 
as it integrates over portions of the plasma with different 
ionization ages. Borkowski, Lyerly & Reynolds (2001) 
showed that NEI models with a single value of r are poor 
descriptions of the so-called "Sedov models" (Hamilton, 
Sarazin & Chevalier 1983), designed to model thermal 
X-ray emission from SNRs. Instead, plane-parallel mod- 
els with multiple temperatures and r are good fits to the 
Sedov models, further justifying the preference for the 
VPSHDCK+VPSHOCK model. We have utilized version 2.0 
of the NEI models in XSPEC, as used by Z06; the files, 
with updated "inner shell processes" (Z06), were kindly 
provided to us by the author of the models, K. Borkowski 
(2007, private communication). 

The data and model fits are plotted in Figs. [3] and 
|U Details of the fits to the EPIC-pn data as well as 
the derived parameters are given in Table [2] For the 
VPSHDCK+VPSHOCK model, the range of ionization ages, 
r = n e ti on , is bounded by < r < t u . Consistent with 
previous studies of SNR 1987A in the X-ray, the upper 
limit to the ionization age, r u , is higher for the lower 
temperature component. The computed abundances are 
given relative to their solar values (Anders & Grevesse 
[1989], hereafter AG89; Wilms, Allen & McCray [2000], 
hereafter WOO; see g]), as defined in XSPEC. Galactic 
foreground absorption is fixed at 6 x 10 20 cm~ 2 , follow- 
ing H06. To account for absorption by the LMC, we take 
the elemental abundances to be 0.5 relative to their solar 
values, with the exception of helium which is kept at its 
solar value. To get a handle on the systematic errors, we 
consider both the AG89 and WOO tables. For meaningful 
comparisons, output abundance values must be normal- 
ized to values relative to hydrogen using the respective 
abundance table used. In the said version of XSPEC, we 
note that the help file for the VPSH0CK function erro- 
neously lists it as being hardwired to the AG89 tables, 
even though it is superseded by the abund command (K. 
Arnaud 2007, private communication). 

The errors are computed using the error function for 
A\ 2 = 2.706. In general, the dependence of the con- 
fidence level on A% 2 is a non-trivial function (Wall & 
Jenkins 2003). Avni (1976) finds that it depends on 
the number of parameters that are estimated simulta- 
neously, i.e., "interesting parameters", and not on the 
total number of parameters in the fitting function. For 
one interesting parameter, the 90% confidence interval 
corresponds to A\ 2 = 2.71; if more interesting parame- 
ters are considered, the required value of A\ 2 increases. 
The Ax 2 = 2.706 level is often quoted as the "90%" con- 
fidence interval in XSPEC analyses of SNR 1987A — we 
wish to point out that this should be regarded with some 
care. Following Z06 and H06, the abundances of He, C, 
Ar, Ca and Ni are held fixed — He and C values are taken 
from Lundqvist & Fransson (1996, hereafter LF96), while 
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Ar, Ca and Ni ones are from Russell & Dopita (1992). 

The VPSHOCK+VPSHOCK model is marginally superior to 
the VNEI +VRAYM0ND one for fitting the EPIC-pn data - 
the reduced chi-square is Xr = 1-21 versus 1.31. For ref- 
erence, we note that the calibration spectra for XMM 
are typically fitted with a "goodness" of x 2 ~ 1-5. To 
check for consistency between the different XMM instru- 
ments, we have taken the VPSHOCK+VPSHOCK fit to the 
EPIC-pn data, folded it with the corresponding detector 
response and compared it with the EPIC-MOS and RGS 
spectra. All of the fit parameters are held fixed in the 
comparison, but the normalization is allowed to vary via 
a constant factor. This yields xl = 1-52 (MOS1), 1.41 
(MOS2) and 1.56 (RGS); for the VNEI+VRAYMOND model, 
we get Xr — 1-45, 1.46 and 1.54 instead. Full-fledged 
and independent fits — like for the EPIC-pn spectra — 
to the EPIC-MOS data and their corresponding Xr val- 
ues are displayed in Fig. |4j We note that making more 
parameters free did not necessarily improve the fits, im- 
plying good cross calibration between the EPIC-pn and 
EPIC-MOS cameras. 

The better counting statistics of the EPIC-pn data 

- and its wider energy coverage compared to the RGS 

- motivate us to use it as a template for obtain- 
ing the temperatures of the two plasma components. 
Moreover, the higher-temperature component lies be- 
yond the energy range of the RGS. From performing 
the VPSHOCK+VPSHOCK fit to the EPIC-pn data, we ob- 
tain fcTi ow = 0.4 keV and /cXhi g h = 3.0 keV, consistent 
with the values of - 0.5 and - 3 keV by Z06 (2004 
Aug/Sept), who also used a two-shock, NEI model. By 
contrast, the VNEI+VRAYMOND fit to the EPIC-pn data 
gives fcT low = 0.3 keV and fcT high = 2.4 keV. An NEI 
plasma is at a higher temperature than is reflected by 
the ionization stages of its elements; therefore, modeling 
it as a plasma in collisional ionization equilibrium under- 
estimates its temperature. We do not expect the temper- 
atures of the plasmas to be substantially lower relative 
to earlier epochs, as the hot spots on the equatorial ring 
continue to brighten. A worthwhile future endeavor will 
be to obtain kT\ ow and fc Thigh for all of the existing Chan- 
dra, XMM and Suzaku observations of SNR 1987A, ana- 
lyzed self-consistently using the same VPSHOCK+VPSHOCK 
model with both AG89 and WOO abundance tables. 

3.2. INTEGRATED FLUX 

By integrating the X-ray spectra, one can obtain fluxes 
in different sub-bands. Park et al. (2004) trisected the 
spectral range into: the 0.3 — 0.8 keV sub-band to rep- 
resent the O line features ("O band"); the 0.8 — 1.2 keV 
sub-band for the Ne line features ("Ne band"); and the 
1.2 — 8 keV sub-band for the Mg/Si lines plus any hard- 
tail emission features ( "H band" ) . The soft X-ray band is 
commonly regarded to be in the 0.5 to 2 keV range (P05; 
H06), a historical convention from the era of ROSAT 
(Hasinger, Aschenbach & Triimper 1996). The "hard 
band" is taken to be 3 to 10 keV by P05. We integrate 
the uncorrected X-ray spectra for various sub-bands and 
tabulate them (for l-cr confidence intervals) in Table [H 
To correct for Galactic and LMC absorption, one sets the 
column densities to zero and then re-computes the fluxes. 
Due to the different abundance tables used and column 
densities inferred, the absorption-corrected luminosities 
are model-dependent. Assuming a distance to the SNR 



1987A of d w 50 kpc, we obtain for L .5-2kcv: 1-95 x 10 36 
erg s" 1 (VPSHOCK+VPSHOCK, WOO tabie), 2.32 x 10 36 erg 
s _1 (VPSHOCK+VPSHOCK, AG89 table) and 2.28 x 10 36 erg 
s" 1 (VNEI+VRAYMOND, WOO table). For L .5-iokcV, we 
get: 2.17 x 10 36 erg s" 1 (VPSHOCK+VPSHOCK, WOO table), 
2.54 x 10 36 erg s" 1 (VPSHOCK+VPSHOCK, AG89 table) and 
2.49 x 10 36 erg s _1 (VNEI+VRAYMOND, WOO table). For the 
2007 Jan observation of P07, i .5-iokoV ~ 2.35 x 10 36 
erg s . 

3.3. LINE FITTING 

For the RGS observations, we consider only data points 
between 0.45 and 1.1 keV, the range over which 11 emis- 
sion lines were previously listed by H06 (see their Table 
5 and also Z06). We note that 38 lines were considered 
by H06, but only 11 were listed so as to facilitate com- 
parison with the Chandra studies. We rebin the spectra 
such that each bin has a minimum of 30 photon counts; 
we find that if the minimum count is 40, some of the 
lines are covered by only 2 to 3 bins. The RGS spectra 
are modeled using a thermal bremsstrahlung component 
for the continuum and a set of Gaussian profiles for the 
lines. We first fit only to the 11 lines of N VII, O VII, 
O vin, Ne ix and Ne x (Fig. HJ. We see that several 
lines are not fitted by such a model — we propose that 6 
additional lines are now observed: the 725, 727, 739, 812 
and 826 eV lines of Fe XVII; and possibly the 873 eV line 
of Fe xviii. (See Table 1 of Behar et al. [2001] and Fig. 
2 of Rasmussen et al. [2001] for details on the expected 
lines in SNR spectra within the RGS energy range.) The 
detection of the iron lines lends support to our derived 
Fe abundance and its implications (see 2]). 

Following H06, we force the widths of all of the lines 
to be the same. The more energetic lines are proba- 
bly powered by the faster shocks, and should therefore 
have larger widths. A more realistic approach is therefore 
to model the widths as some increasing function of the 
ionization potential of each line. However, the rather 
narrow energy range considered here does not warrant 
such an approach. The redshift inferred from the fits is 
(1.88 ± 0.26) x 10" 3 . The Gaussian width of the lines is 
erfit = 0.87 ±0.11 eV — the full width at half-maximum 
(FWHM) value is 2.05 ±0.26 eV, which is narrower than 
the 5.3 ± 1.0 eV value of H06. The measured Gaussian 
width is indicative of the velocities of the shocks respon- 
sible for the X-ray emission — the velocities inferred are 
v s ~ 0.87c/ E c y, where E e y is the energy of a given line 
in electron volts. For a 1 keV line, we have v s ~ 300 
km s _1 , consistent with the Chandra studies of Z06, who 
find that the X-ray emission originate from shocks with 
300 < v s < 1700 km s -1 . 

The fluxes we obtain from fitting to these lines are 
given in Table O In Fig. [SJ we show the relative in- 
creases in flux of the lines (or line complexes) measured 
in our study versus those of H06. Since 2003 May, the 
line fluxes have increased by factors ~ 3 to 10 — the 775 
eV line of O vin shows the strongest increase. On aver- 
age, the lines show an increase of 6.0 ±0.6 over ~ 4 years. 
The 2003 May soft X-ray flux was measured by H06 to 
be ^ .5-2kcV = (8.10 ±0.09) x 10~ 13 erg cm" 2 s _1 . Our 
2007 Jan value is i^ .5-2kcV = (3.34 ± 0.04) x 10~ 12 erg 
cm~ 2 s _1 , corresponding to an increase of 4.1 ±0.1 in 
the soft X-ray flux. Generally, the same shocks that are 
responsible for the brightening of the optical and soft X- 



4 



Heng et al. 



ray hot spots are also powering the lines. The differences 
in the flux increases may be an indication that the equiv- 
alent widths of the lines are changing, which is evidence 
for the evolution of the elemental abundances. Such an 
investigation is deferred to future studies, when one can 
obtain smaller error bars for the relative increases of the 
line fluxes (Fig. [S]). 

4. DISCUSSION 

4.1. INDIVIDUAL ABUNDANCES & THEIR RATIOS 

We compare the derived elemental abundances to those 
of Z06 and H06 and list them relative to hydrogen in 
Table [3l uncertainties in the AG89 and WOO abundance 
tables are not propagated. We emphasize that care must 
be taken to specify the abundance table used, as this may 
lead to widely differing values of the derived abundances 
(relative to hydrogen). In modeling the LMC absorption, 
Z06 used the elemental abundance table of AG89, while 
H06 chose the table of WOO because the lower oxygen 
abundance fitted the K absorption edge in the EPIC data 
better. 

Fransson et al. (1989) found N/O = 1.6 ± 0.8, about 12 
times higher 3 than the AG89 solar value, which they in- 
terpret as evidence of substantial CNO processing. LF96 
found N/O = 1.1 ± 0.4, while Sonneborn et al. (1997) 
found N/O = 1.7 ± 0.5. All three of these values were 
derived from optical/ultraviolet data. Next, we turn our 
attention to the N/O ratio derived from X-ray stud- 
ies. Linearly propagating the errors listed by Z06, we 
find that their results yield N/O = 1.10±°;45; tne y re- 
mark that their derived C, N and O abundance, C+N+O 
w 1.98 x 10~ 4 , is lower by about a factor of 2 compared to 
the 3.72 xl0~ 4 value of LF96. Model A (VNEI+VRAYMOND) 
in H06 yields C+N+O w 1.67 x 10~ 4 and a rather wide 
range in the nitrogen-to-oxygen ratio, N/O = 1.33lJ'^. 
Our VPSHOCK+VPSHOCK fit to the EPIC-pn data yields 
C+N+O re 1.29 x 10~ 4 and N/O = IA7±° 3 3 1 (using the 
WOO table); we call this combination of N and O the 
"best fit point" . Note that the C abundance is held fixed 
at 0.09 relative to solar (~ 3 x 10~ 5 relative to hydrogen) 
for the C+N+O values derived from the X-ray studies. 

We next perform a more careful analysis of the N/O 
ratio. We first generate a \ 2 map quantifying the inter- 
dependence of the fits to the N and O abundance. Con- 
tour lines in the \ 2 map form "error ellipses" , which are 
shown for different A% 2 values from the best fit point 
(Fig. ED. At the Ax 2 = 2.706 level, N/O = 1.17 ± 0.20 
for the EPIC-pn data. In linearly propagating the errors 
in the individual abundances, one is in essence adopting 
the largest possible range of ratios, which can be visu- 
alized as the edges of a rectangle in the contour map. 
Our error analysis improves the uncertainties because it 
considers only values of the abundance ratio within the 
specified contour. Within the same confidence interval 
considered, the corresponding C+N+O value is from ~ 
1.1 to 1.4 x 10~ 4 . We see that the errors in the N/O 
ratio and the individual N and O abundances can be 
constrained at the ~ 20% level. We thus confirm the 
C+N+O under-abundance noted by the Chandra studies 
of Z06, who suggest a couple of physical reasons for such 
a result: the sub-LMC abundance of C+N+O within 

3 The AG89 and WOO values for solar nitrogen-to-oxygen abun- 
dance are N/O = 0.132 and 0.110 by number, respectively. 



the progenitor star, Sanduleak -69°202, and/or an extra 
source of possibly non-thermal, X-ray continuum. 

We perform the same analysis for N/S (Fig. [5D. Again 
using the WOO table, linear propagation of the errors in 
the abundances obtained from the VPSHOCK+VPSHOCK fit 
yields N/S = 4.59t 1 'gg, while the error ellipse analysis 
gives N/S = 4.59^ 1 '33- The iron, nitrogen and oxygen 
lines are predominantly from the lower energy part (< 1 
keV) of the spectrum, while the sulphur lines are situ- 
ated between ~ 2.2 and 3.1 keV. (The silicon lines are 
located between ~ 1.8 and 2.2 keV.) Abundance ratios 
based on lines of widely differing energies lead to rounder 
error ellipses — "error circles" . In such cases, a more 
thorough error analysis will not constrain the abundance 
ratio better, as in the case of N/S. This is partially an 
instrumental effect — when the energy resolution is com- 
parable to the spacing of the lines, the line complexes 
overlap and are only partially resolved. There is also the 
issue of choosing a coordinate system in which the fit- 
ting parameters are "orthogonal" . Hydrogen and helium 
have no X-ray lines — their abundances are derived from 
the strength of the X-ray continuum. When the pair of 
elements considered are situated close to each other in 
energy, the hydrogen abundance has to first order a lin- 
ear dependence on the continuum strength. Thus, the 
ratio of the considered abundances is tightly constrained 
as the individual abundances track each other closely. 
By contrast, when the pair of abundances considered are 
located far apart in energy, this linear dependence of hy- 
drogen abundance on the continuum is broken as the 
dominant uncertainties are in temperature rather than 
in flux. Orthogonality is now absent. To attain orthogo- 
nality for such pairs of lines, one has to construct models 
that directly fit to the abundance ratio considered, an 
approach which is not explored in this paper. 

The error in the abundance ratios increases as one 
moves a given A\ 2 away from the minimum point. In 
Fig. [9l we compute the mean error sustained by the 
various ratios as a function of A% 2 . According to Avni 
(1976), if three interesting parameters are considered 
(see E)3.ip . the 90% confidence interval is situated at 
Ax 2 = 6.25. In this case, the N/O abundance ratio 
suffers from errors ~ 25%. 

4.2. THE PROGENITOR OF SNR 1981 A: 
SINGLE-STAR OR BINARY MODEL? 

A more revealing approach to analyze the elemental 
abundances is to normalize the results listed in Table[3]by 
the "canonical" values of the LMC abundances (Hughes, 
Hayashi & Koyama 1998). These were derived using a 
sample of 7 middle-aged SNRs in the LMC (N23, N49, 
N63A, DEM 71, N132D, 0453-68.5 and N49B). This ap- 
proach was first explored by A07, who showed that the 
normalized abundances appear to cluster in two groups: 
N, O, Ne, Mg and Fe are slightly more than half their 
LMC values, while Si, S and Ni exceed their LMC values. 

We generalize the A07 approach by considering both 
sets of abundances derived from using the AG89 and WOO 
tables. The normalized abundance, relative to its respec- 
tive LMC value, is i?87A/LMci it is plotted as a function 
of the elemental mass number in Fig. 1101 The error 
bars for i?87A/LMC are computed by linearly propagating 
the errors listed in Table [3] and in Hughes, Hayashi & 
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Koyama (1998). We caution that additional systematic 
errors may be present that are not taken into account. 
For example, the derived abundance for Fe is a sensi- 
tive function of temperature and may vary substantially 
when small changes are made to Xi ow . 

We see that the elements O, Ne, Mg and Fe are under- 
abundant, while Si and S are over-abundant, consistent 
with the findings of A07. With the exception of Fe, there 
is a tendency for i?s7A/LMC to increase with larger ele- 
mental mass number, a trend that is independent of the 
AG89 or WOO tables, though we note that it is more pro- 
nounced with the latter. The Fe abundance derived is 
essentially independent of the AG89 or WOO tables, and 
is under-abundant by about 70% relative to the LMC 5 . 
The under-abundance of Fe and O was previously noted 
by Hasinger et al. (2006), who argued for the existence of 
iron-oxygen "rust grains" . The reduced abundance of Fe 
alone suggests that the iron is locked up in dust grains. 
However, Dwek & Arendt (2007) showed from an analy- 
sis of the infrared-to- X-ray flux ratio — TZmx < 1 versus 
the theoretically expected value of ~ 10 2 to 10 3 — that 
the dust in SNR 1987A is severely depleted compared to 
standard dust-to-gas mass ratios in the LMC, suggesting 
low dust condensation efficiency or dust destruction in 
the hot X-ray gas. In fact, TZmx was shown to decrease 
with time, which is direct evidence for dust destruction. 
Our derived plasma temperatures are consistent with this 
scenario — even if dust could form, it would be destroyed 
at these temperatures. 

In light of Fig. [T0l the central question to ask is 
whether the progenitor of SNR 1987A arose from a single 
star or a binary system? Sanduleak -69°202 was known 
to be a blue supergiant (BSG) at the time of the super- 
nova explosion, contrary to the expectation that massive 
stars end their lives as red supergiants (RSGs). Observa- 
tions of low- velocity, nitrogen-rich circumstellar material 
are interpreted as the progenitor star being a RSG until 
about ~ 20, 000 years before its death (Fransson et al. 
1989). Such a time scale has in turn been interpreted as 
the Kelvin-Helmholtz time of the helium core (Woosley 
et al. 1997). This BSG-RSG-BSG evolution remains 
one of the greatest challenges for the single-star model 
(Woosley et al. 1997; Woosley, Heger & Weaver 2002), as 
is the observed system of three rings produced ~ 20, 000 
years before the explosion. The favored single-star mod- 
els require a combination of reduced metallicity and "re- 
stricted semi-convection" (Woosley 1988), the former of 
which is supported by our derived abundance ratios. An 
additional challenge for the single-star model is to repro- 
duce the over-abundance of Si and S (Fig. [10|), which 
may require the invoking of some non-standard mixing 
process (H.-T. Janka 2007, private communication). 

Binary solutions to Sanduleak -69° 202 are sub-divided 
into accretion and merger models (see Woosley, Heger & 
Weaver [2002] and references therein). The binary ac- 
cretion models allow helium- and nitrogen-rich material 
to be added to the progenitor star and require the disap- 
pearance of the mass donor in an earlier supernova event. 
In the binary merger scenario proposed by MP0607, two 

4 The Fe abundance obtained from the VPSHOCK+VPSHOCK fit 
(WOO table) varies by ~ 16% when Ti ow is changed by ~ 10%. 

5 The O abundance relative to H for the AG89 and WOO tables 
are the same. 



stars with masses ~ 5Mq and ~ 15Mq are initially or- 
biting each other with a period ~ 10 years. The more 
massive companion transfers mass to the less massive 
star only after the former has completed helium burn- 
ing in the core. A common envelope (Paczyhski 1976) is 
formed, during which core material from the primary star 
is dredged up to the surface. The merger process takes a 
few hundred years, culminating in an initially over-sized 
RSG, which loses its excess thermal energy over a few 
thousand years to become a BSG. The spun-up, rapidly 
rotating BSG produces a fast stellar wind, sweeping up 
ejecta associated with the merger, producing the triple 
ring nebula we now see in projection. The nearly axi- 
symmetric but highly non-spherical nature of the rings 
suggests that rotation played a role in their formation 
and is consistent with the proposed scenario. The beauty 
of the model lies in the fact that it requires no physically 
ad hoc assumptions — apart from a small kick of ^ 2 km 
s" 1 given to the ejecta to displace the center of the outer 
rings from their symmetry axis — and makes a number of 
predictions. In their favored model, MP0607 assert that 
the outer rings are ejected before the stellar core material 
is dredged up, while the inner, equatorial ring — the site 
of the observed X-ray emission — is ejected afterwards 6 . 
A relevant consequence of this model is that the dredged- 
up heavy elements will manifest themselves in the form 
of X-ray emission lines. This may explain the trend we 
see in Fig. [TU] — the challenge for binary merger models 
is to reproduce the derived i?87A/LMC values. 
Stellar nucleosynthesis can add and not subtract iron 
in the context of the MP0607 binary merger model, 
we expect i?87A/LMc(F e ) > 1- We are again led to the 
question: where is the iron? If the Fe abundance derived 
is an upper limit on the iron used to form Sanduleak 
-69°202, then it is clearly inconsistent with "standard" 
Fe abundances in the LMC. An alternative interpreta- 
tion is that there is a strong spatial variation in the Fc 
abundance throughout the LMC, such that the iron is 
sub-LMC at the site of SN 1987A and equal to its LMC 
abundance elsewhere. The C+N+O under-abundance 
suggested in §4.11 supports such a conclusion. An im- 
provement over using the Hughes, Hayashi & Koyama 
(1998) ASCA abundance values is to re-analyze and ex- 
pand upon their SNR sample using Chandra and XMM. 
Existing studies (e.g., Hughes et al. 2006) tend to pick 
out regions of interest that may include ejecta enrich- 
ment; instead, X-ray emission should be extracted from 
the entire blast wave, from which average abundances 
can be inferred. Future studies will be invaluable to- 
wards resolving these issues. 
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6 As noted by MP0607, this hypothesis is verifiable/refutable, as 
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FlG. 1. — Schematic representation of the physical configuration in SNR 1987A. Courtesy of Richard McCray (JILA, University of 
Colorado) and the Chandra press release team. Figure has been omitted due to large file size; please refer to the electronic version of the 
Astrophysical Journal for this figure. 
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TABLE 1 

XMM-Newton observations of SNR 1987A 



Instrument 


Read-out 


Filter 


Satellite 


Date 




Time (UT) 


Net Exposure 




Mode 




Revolution 






(s) 


EPIC-pn 


FF, 73 ms 


Medium 


1302 


2007 Jan 17- 


-19 


19:28—01:18 


93723 


EPIC-MOS1 


FF, 2.6 s 


Medium 


1302 


2007 Jan 17- 


-19 


18:24—01:18 


109479 


EPIC-MOS2 


FF, 2.6 s 


Medium 


1302 


2007 Jan 17- 


-19 


18:24—01:18 


109560 


RGS1 


Spectro 




1302 


2007 Jan 17- 


19 


18:23—01:19 


110528 


RGS2 


Spcctro 




1302 


2007 Jan 17- 


-19 


18:23—01:19 


110475 



TABLE 2 

Comparison of the VNEI+VRAYM0ND and the VPSH0CK+VPSH0CK models 



Parameter 


VNEI+VRAYM0ND 


VPSH0CK+VPSH0CK 


VPSH0CK+VPSH0CK 




WOO* 


AG89* 


WOO* 


x 2 


820.2032 


805.6812 


755.9608 


xt 


624 


623 


623 


1.314428 


1.293228 


1.213420 


<w ( k °V) 


S20+ 010 
u.ozu_ 00g 


0.419lH 2 6 o 


434+0-018 
u - 4 ' r ' -0.008 


^high (keV) 


9 49 +0.175 
-0.176 




2 99+ ' 24 
^• aa -0 .22 


Galactic N H (10 21 cm- 2 )t 


0.6 


0.6 


0.6 


LMC Nn (10 21 cm- 2 ) 


3.13 ± 0.14 


2.37 ±0.10 


2 56+ ' 07 
Z - OD -0.11 


Tu(T low ) (10 11 cm" 3 s) 




o fi ,+1.60 


O OO + 1.02 


Tu(T Mgh ) (10 11 cm- 3 s) 


1 in+0.31^ 
1 - ±a -0.24 


2 46+°' 79 


2 21+ 1 ' 44 


Het 


2.57 


2.57 


2.57 


ct 


0.09 


0.09 


0.09 


N 


116+ 0172 

U.11D_q 1()6 


0.303 ± 0.060 


^fi+0.108 


O 


0228+ ' 0074 
u.uzzo_ 0058 


0475+0-0045 
u.uko_ 0044 


0502+ (UIH:>s 
u.uouz_ 0060 


Ne 


n i c 7 +o.oi7 

U - i0 '-0.015 


156+ ' 012 


204+ ' 017 
u.zu*_ 01g 


Mg 


208+ 021 


116+ 007 

U.11D_ 00g 


168+ ' 011 


Si 


0-743±g : g|| 


262+ ' 025 
u.zdz_ 017 


479+0-046 
—0.030 


s 


0.689l° o:il ? 


452+ ' 059 


^85+ ' 075 
u.ooo_ 074 


Art 


0.54 


0.54 


0.54 


Cat 


0.34 


0.34 


0.34 


Fe 


n-io.7+0-0038 
U.Ujy ( _ 0033 


0630+ 0041 

U.UDJU_g 001 g 


n n8q5+ '00 50 

u.uoyo_ 0060 


Nit 


0.62 


0.62 


0.62 



Note: Abundances of elements arc listed relative to solar, 
f Parameter is held fixed (see text). 

t. In the VNEI+VRAYM0ND model, the low and high temperature components arc from the VRAYM0ND and VNEI fits, respectively. 
£ Number of degrees of freedom in fit. 

Abundance table in XSPEC from Anders & Grcvesse (1989). 

* Abundance table in XSPEC from Wilms, Allen & McCray (2000). 

\? VNEI model considers only one value of r — r u . 



Fig. 2. — XMM-Newton EPIC image of the region around SNR 1987A. The red-green-blue (RGB) color image is composed of images 
from three energy bands (0.2 — 1.0, 1.0 — 2.0 and 2.0 — 4.5 keV) and from all three of the EPIC instruments (PN, MOS1 and MOS2). The 
individual images arc exposure-corrected and out-of-time event-subtracted (for EPIC-pn). The source and background extraction regions 
used for the spectral analysis are denoted by "1987A" and "bg", respectively. The blue, arc-like structures are due to photons from the 
bright X-ray binary LMC X-l (located outside the field-of- view) , which are singly-reflected on the mirror shells. Figure has been omitted 
due to large file size; please refer to the electronic version of the Astrophysical Journal for this figure. 
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TABLE 3 

Elemental abundances from different studies (relative to hydrogen) 



Element 



Z06^ 
VPSHOCK+VPSHOCK 



H06* 

VNEI+VRAYMOND 



This Study <•> 
VPSHOCK+VPSHOCK 



This Study* 
VPSHOCK+VPSHOCK 



(7.83ti;S) 

(3.57tg;S)_ 
(9.12 ± 1.52). 



s-oo±i:»l) 

4-27±8;S) 
2.51^) 



N 

O 

Ne 

Mg 

Si 

S 

Fe 



9.93 



7.30Z 



+1.42 > 



(7.48t°;3l) 



(3.40 ± 0.67) _ 




o 95 +o.i9> 



6.54 



-0.43 ' 



1.70] 



(1.09±0.14)_ B 

(2.83t°;J 6 9 )_ 6 



Note: A& is shorthand notation for A X 10 . 

Abundance table in XSPEC from Anders & Grevesse (1989). 

4k Abundance table in XSPEC from Wilms, Allen & McCray (2000). 



TABLE 4 

AbSORPTION-UNCORRECTED EPIC-pn FLUXES IN VARIOUS SUB-BANDS 

Sub-band photons cm~ 2 s _1 erg cm~ 2 s _1 
0.2 — 0.8 kcV (1.16 EE 0.01)_ 3 (1.14 EE 0.01)_ 12 

11+0.02 \ ' (, c -,+0.03> 



0.8 — 1.2 kcV 



-3 



-12 



1.2 — 8 kcV (4-93l°;J 5 

0.5 — 2 kcV (2.28 ± 0.02)_ 3 (3.34 ± 0.04) 

3-10keV (5.3ll° :S)_ B (3.83l°; 23 ) 

0.5 — 10 kcV (2.42+° ^) (4.05 ± 0.07) 

fo fi-7+O-OsS 



-13 

-12 



0.2 — 10 kcV 



-3 



(4.24±g;gg)_ 



Note: 1-(T (68%) confidence intervals; Ab is shorthand notation for A X 10 . 



TABLE 5 
Line fluxes from RGS observations 



Line 


Energy 


Flux 




(eV) 


(10~ 5 photons cm~ 2 s _1 ) 


N vn 


500 


18.67ti;S 


O VII 


561 


4.94+i^ 


O VII 


569 


1 o ,+ 2 - 12 

1 - uo -1.05 


O VII 


574 


8-92_ 2 ; 37 


O VIII 


654 


25.18± 2 /^ 


O VII 


666 


q ill +0.65 

a.zi_ 1 37 


Fe xvii 


725 


00+0.74 

3-83_ 1 g3 


Fe xvii 


727 


12.98t 2 ;J 5 


Fe xvii 


739 


7 19+0. 63 


O VIII 


775 


9.40t;-«| 


Fe xvii 


812 




Fe xvii 


826 


16.001^3 


Fe xviii 


873 


6-00±S;IS 


Ne ix 


905 




Ne ix 


915 


1 q 1 +2.58 


Ne ix 


922 


18.87t 2 ;gg 


Ne x 


1022 


17.75+^f 


O VII 


561+569+574 


14 -9 1 + 5'l5 


Ne ix 


905+915+922 


29.64+fsi 



Note: Ax 2 — 2.706 confidence intervals. 
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2007 XMM-NEWTON EPIC-pn: SNR 1987A 
1 2 1 ' ' — ' — ' — ' ' ' 1 1 ' ' — ^~ 




0.3 0.4 0.5 0.6 0.7 0.80.9 1 2 3 4 567 910 



Channel Energy (keV) 



Fig. 3.— XMM EPIC-pn data, fitted by the VNEI+VRAYMOND and VPSHQCK+VPSHQCK models in XSPEC. The insert zooms in on a possible 
detection of the Fe K« line. Details of the model fits, as well as the derived parameter values, are given in Table [2] 
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2007 XMM-NEWTON EPIC-MOS1: SNR 1987A 



10' 



1 10° 



10" 



I ^o- A 

o 



10"° 
6 
4 
2 

x 
-2 
-4 
-6 




Data 

VPSHOCK+VPSHOCK fit 
VNEI+VRAYMOND fit 

— , , , — , — , — | 



if* 





0.3 0.4 0.5 0.6 0.7 0.80.9 1 2 3 4 5 6 7 

Channel Energy (keV) 




Fig. 4.— XMM EPIC-MOS1 (top) and EPIC-MOS2 (bottom) data, fitted by the VPSHQCK+VPSHDCK model. For the VPSHOCK+VPSHOCK 
model, the reduced chi-square values are \r = 1-53 (MOS1) and 1.40 (MOS2). For the VNEI+VRAYMOND model, the reduced chi-square values 
are \l = 1-45 (MOS1) and 1.46 (MOS2). 



XMM-NEWTON 87A 



11 



2.0 



1.5 



VIII 



RGS 1 
RGS 2 




/, ' Wl 



0.7 0.8 0.9 

Channel Energy (keV) 



1.1 



2.0 



1.5 



1.0 



E 0.5 
o 



VIII 



RGS 1 
RGS 2 



N VII 



Fe XVII Fe XVII 



VII 



VII 



VIII 



Fe XVIII 



Ne IX 

| Ne X 



WW 




0.7 0.8 0.9 

Channel Energy (keV) 



1.0 



1.1 



Fig. 5.— XMM RGS data considering 11 (top; \l = 3-37 for JV do f = 652) and 17 (bottom; Xr = 1 - 92 for ^dof = 645) Gaussian fits to 
the lines, respectively (see text). 
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E 6 - 



500 600 700 800 

Line energy (eV) 



900 



1000 



Fig. 6. — Flux increase in emission lines from XMM RGS data between 2003 May (Haberl et al. 2006) and this study (2007 Jan). The line 
complexes (denoted "LC") represent the summed fluxes from three lines. See Tableland the text for details of the lines. The horizontal, 
dotted line shows a fit for the average increase in the line fluxes: 6.0 ± 0.6 (x, 2 = 0.98). 




0.03 0.04 0.05 0.06 0.07 

abundance (relative to solar) 



Fig. 7. — Contour map of N and O for various Ax 2 values from the best fit point, marked by the large cross. Values of N and O residing 
along the dotted line have the same N/O value as that for the best fit point. 
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Fig. 8.— Same as Fig. [7] but for N/S. 
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Anders & Grevesse (1989) O 
Wilms, Allen & McCray (2000) A 



1.5 - 



1.0 -■ 



Ne 



Si 
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<> 
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Fig. 10. — Abundance ratios of elements found by our study of SNR 1987A to the LMC quantities derived by Hughes, Hayashi & Koyama 
(1998). In our XSPEC analysis, we examine abundances obtained using the tables of Anders & Grevesse (1989) and Wilms, Allen & McCray 
(2000). 



